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System of two quantum dots attached to external electrodes is considered theoretically in orbital 
Kondo regime. In general, the double dot system is coupled via both Coulomb interaction and 
direct hoping. Moreover, the indirect hopping processes between the dots (through the leads) are 
also taken into account. To investigate system's electronic properties we apply slave-boson mean 
field (SBMF) technique. With help of the SBMF approach the local density of states for both 
dots and the transmission (as well as linear and differencial conductance) is calculated. We show 
that Dicke- and Fano-like line shape may emerge in transport characteristics of the double dot 
system. Moreover, we observed that these modified Kondo resonances are very susceptible to the 
change of the indirect coupling's strength. We have also shown that the Kondo temperature become 
suppressed with increasing asymmetry in the dot-lead couplings when there is no indirect coupling. 
Moreover, when the indirect coupling is turned on the Kondo temperature becomes suppressed. 
By allowing a relative sign of the nondiagonal elements of the coupling matrix with left and right 
electrode, we extend our investigations become more generic. Finally, we have also included the level 
renormalization effects due to indirect tunneling, which in most papers is not taken into account. 
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I. INTRODUCTION 

Originally, Kondo effect was discovered in non- 
magnetic metal containing magnetic impurities at low 
temperature. The effect comes from strong electron cor- 
relations and can be regarded as interactions of the impu- 
rity spin with cloud of the conduction electrons in metal. 
Scattering of the conduction electrons from the localized 
magnetic impurities leads to increase of the resistivity at 
low temperature X In recent two decades one could no- 
tice revival of this effect as it was predicted and observed 
in transport through quantum dotsi^r— Although, in this 
case Kondo resonance leads to increasing conductivity 
(not resistivity as in original Kondo effect) with decreas- 
ing temperature below so-called Kondo temperature Tk, 
the physical mechanism of the phenomena is common. 
Here, the role of the magnetic impurity plays a spin on 
the dot. 

However, despite of regarding spin degree of freedom 
one may assume any two- valued quantum numbers, as for 
instance, an orbital degree of freedom, to realize Kondo 
effect. In the case of (at least) two discrete orbital levels 
coupled to the external leads one may deal with orbital 
Kondo effect To explain mechanism of creation of the 
Kondo state let us introduce spinless electrons in a sys- 
tem of two single-level quantum dots coupled to external 
leads as shown in Fig. [TJ The dots' energy levels are well 
below the Fermi level of the leads. Due to large inter- 
dot Coulomb repulsion only one of them is occupied by 
an electron. On the other hand, removing an existing 
electron on the double dot requires adding energy to the 
system. Thus, the system is in the deep Coulomb block- 
ade and the sequential tunneling events are prohibited. 
However, according to the Hcisenberg uncertainty princi- 
ple, the higher-order processes, however, may appear on 



a very short time scale. Assuming that initially the upper 
dot was occupied, then it can tunnel onto the Fermi level 
of the lead (left or right) and simultaneously another elec- 
tron from the Fermi level of the left or right electrode may 
tunnel to the lower dot. As a result, the charge exchange 
occured between the dots. A coherent superpostion of 
such coherent events give rise to the sharp resonance in 
the density of states at the Fermi level. 

Altough orbital Kondo effect has been investigated 
in different geometries of the two orbital level system 
both experimentally^ and theoreticallj*i£~— there is a 
little comprehensive studies on the influence of the indi- 
rect coupling and/or asymmetry couplings on this phe- 
nomenon. In the case of spin Kondo effect most of re- 
searchers assume the maximal value of this coupling 
However, the maximal value of the indirect coupling in 
the case of the orbital Kondo effect may leads to suppres- 
sion of Kondo resonance. Specifically, this effect is de- 
stroyed when the magnetic flux achieves 2n7r (n £ Z) 23 i 24 
due to the formation of the bound state in the continuum 
(BIC)2£. Recently, Kubo et al. have investigated both 
spin and orbital Kondo effect in double quantum dot re- 
garding non-maximal values of the indirect couplin g i 26 . 
They have found that in the condition of intermediate in- 
direct coupling strength the differential conductance re- 
veals two kinds of peaks. Altough, we considered similar 
system, the role of this paper is to give more insight into 
the connection between the role of the indirect coupling 
and various interference effects. Moreover, we assume 
different asymmetry coupling of the dots to the external 
leads. More specifically, in the system considered here 
one dot is coupled to the leads with constant strength, 
whereas the coupling of the second dot to the leads can 
be continuously tuned. Time reversal symmetry implies 
that the amplitude of the indirect coupling strength may 
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change signal. Thus, we also include this case in our 
consideration. 

In very recent experiment Tarucha et al^ observed 
that the period of the Aharonov-Bohm oscillations are 
halved and the phase changes by half a period for the 
antibonding state from those of the bonding. They con- 
clude that these features can be related to the indirect 
interdot coupling via the two electrodes. 

We also point out that similar double dot system 
has been investigated in Refs. (29i[30| where the authors 
have found a novel pair of correlation-induced resonances 
as well as they studied the charging of a narrow QD 
level capacitivcly coupled to a broad one. However, the 
correlation-induced effects are not related to the Kondo 
physics. 

As we mentioned before various quantum interference 
effects, which were previously reported in atomic physics 
or quantum optics^Ir— , were also discovered in elec- 
tronic transmission through QDs systems attached to 
the leads^ir— . Here, we consider Dicke and Fano res- 
onances. In original Dicke effect one observes in spon- 
taneous emission spectra a strong and very narrow reso- 
nance which coexist with much broader line. This occurs 
when the distance between the atoms is much smaller 
than the wavelength of the emitted light (by individual 
atom). The former resonance, associated with a state 
which is weakly coupled to the electromagnetic field, is 
called subradiant mode, and the latter, strongly coupled 
to the electromagnetic field refers to superradiant mode. 
In the case of electronic transport in mcsoscopic systems 
(for instance, quantum dots) this effect is due to indirect 
coupling of QDs through the leads. Generally, indirect 
coupling can lead to the formation of the bonding and an- 
tibonding states. As a result, a broad peak corresponding 
to the bonding state and a narrow one referring to the 
antibonding state emerge in the density of state o 39 ' 40 . 

Let us now tell something about Fano effect. In 
experiment, Fano effect manifests itself as asymmetric 
line shape in emission spectra. It comes from quan- 
tum interference of waves resonantly transmitted through 
a discrete level and those transmitted nonresonantly 
through a continuum of states. The effect was ob- 
served in optics and also in electronic transport through 
QDs systems^—. However, in this case the Fano phe- 
nomenon is due to the quantum interference of electron 
waves transmitted coherently through the dot and those 
transmitted directly between the leads 4 ^. 

It is convenient to associate resonant channel with dis- 
crete level and nonresonant channel with continuum of 
states. When the electron wave passes through resonant 
channel its phase changes by 7r (within T), whereas the 
phase of electron waves in nonresonant channel changes 
very slowly around the resonant level (r is the width 
of the discrete level). (Of course, Fano effect occurs 
only when discrete level is embedded into a continuum.) 
Consequently, on the one side of the discrete level elec- 
tron waves through two channels interfere constructively, 
whereas on the other side they interfere destructively. As 




FIG. 1: Schematic picture of the double dot system. The 
parameter a takes into account difference in the coupling of 
the two dots to external leads (a G (0, 1}). Tuning parameter 
a one can change geometry of the system from the parallel 
one for a = 1 to the T-shaped geometry for a = 0. 



a result, one observes asymmetric line in conductance 
around the discrete level position. 

This effect can be also observed in system of two quan- 
tum dots embedded in two arms of AB rin g 21 ' 47 " — or in 
the so-called T geometr y 45 ! 51 . Here, very narrow (broad) 
level, which is weakly (strongly) coupled to the leads, 
corresponds to the resonant (nonresonant) channel. The 
narrow level must appear within the broad one. As men- 
tioned before the difference in the coupling strengths of 
the bonding and antibonding states is due to indirect cou- 
pling. The phase shift of wave function in the broad level 
is negligible when the energy changes within the narrow 
level and Fano resonance may appear. 

In this paper the orbital Kondo effect in electronic 
transport through two coupled quantum dots is consid- 
ered theoretically. Generally, the quantum dots may in- 
teract via both Coulomb repulsion and hopping term. 
To calculate local density of states (LDOS) for both 
dots, transmission, and differential conductance we em- 
ploy slave-boson mean field approach. To show the for- 
mation of the bound state in the continuum as the indi- 
rect coupling strength approaches its maximal value we 
calculate the Fricdel phase. Due to emergence of the BIC 
the Fiedel phase, usually continuous, changes abruptly 
at the energy corresponding to the BIC. Finally, to be 
more familiar with experiment we show differential con- 
ductance. 

The paper is organized as follows. In Section 2 we de- 
scribe the model of a double-dot system which is taken 
under considerations. We also present there the slave- 
boson mean-field technique used to calculate the basic 
transport characteristics. Numerical results on the or- 
bital Kondo problem are shown and discussed in Section 
3. Final conclusions are presented in Section 4. 
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II. THEORETICAL DESCRIPTION 
A. Model 

To investigate various interference effects in Kondo 
regime we consider (spinless) Anderson Hamiltonian for 
double quantum dots coupled to external leads. Exper- 
imentally, such system may be realized by applying a 
magnetic field which lifts spin degeneracy of each dot. 
Generally, this Hamiltonian consists of three parts, 



H = H r + H 



DQD 



(1) 



where the first term, H c , describes nonmagnetic elec- 
trodes in the non-interacting quasi-particle approxima- 
tion, H c = H l + Hr, with Hp = X)k PkpcLpCkp ( for elec_ 
trodes (3 = L,R). Here, (ck/?) creates (annihilates) 
an electron with the wave vector k in the lead j3, whereas 
£k/3 denotes the corresponding single-particle energy. 

The next term of Hamiltonian ([1]) describes two cou- 
pled quantum dots, 



H 



DQD 



where rij = d\di is the particle number operator, Si is the 
discrete energy level of the i-th dot {% = 1,2), t denotes 
the inter-dot hopping parameter (assumed real), whereas 
U is the inter-dot Coulomb integral. 

The last term, Ht, of Hamiltonian ([1} describes elec- 
tron tunneling between the leads and dots, and takes the 
form 



£ l d\d l +t{d\d2 + h.c.) + Un 1 n 2 , (2) 



Hn 



= E E K c U 

W/3 i=l,2 



h.C.), 



(3) 



where V^ k are the relevant matrix elements. 

Finite widths of the discrete dots' energy levels come 
from coupling to the external leads and may be expressed 
in the form T^(e) = 2ir\V.^\ 2 p, where p denotes den- 
sity of states in the left and right lead (px. = pu = p). 
Furthermore, we assume that F^ is constant within the 
electron band, r^(e) = = const for e € (-D, D), and 
rfi( £ ) = otherwise. Here, 2D denotes the electron band 
width. 

For the system taken under considerations, the dot- 
lead couplings can be written in a matrix form, 



r /3 r /3 
1 11 1 12 

22 



1 21 1 



(4) 



where the off-diagonal matrix elements are assumed to 

34^52 Qff.^jagonai ma trix 



be r$ 



-/3 



r F /3 
1 ll 1 22' 



12 _ x 21 — ' 

elements of T 13 take into account various interference ef- 
fects resulting from indirect tunneling processes between 
two quantum dots via the leads. These off-diagonal ma- 
trix elements may be significantly reduced in comparison 



to the diagonal matrix elements or even totally sup- 
pressed due to complete destructive interference. To take 
all those effects into account, the parameters and qn 
are introduced. Furthermore, we assume that qp are real 
numbers and obey the condition \qp\ < 1. When qp is 
nonzero, the processes in which an electron tunnels from 
one dot to the /3-lead and then (coherently) to the an- 
other dot are allowed. Introducing parameter a. which 
takes into account difference in the coupling of the two 
dots to external leads, the coupling matrix Eq. §E§ can 
be rewritten in the form, 



qp\foi 1 



(5) 



Tuning parameter a one can change geometry of the sys- 
tem from the parallel one for a = 1 to the T-shaped ge- 
ometry for a = 0. In the T-shaped geometry upper dot 
is disconnected from the leads. All intermediate values 
of the a refer to an intermediate geometry where each of 
the two dots is coupled to leads with different strength. 
We further assume symmetric coupling T L = r R = T/2, 
where T is energy unit. 



B. Method 

We perform our calculations in large intcrdot charg- 
ing energy limit, more specifically, when U — > oo. 
Slave boson approach is one of the techniques which 
allows investigate strongly correlated fermions in low 
temperatures^. This method relies on introducing aux- 
iliary operators for the dots and replacing of the dots' 
creation and annihilation operators by fjb, fet/j respec- 
tively. Here, slave-boson operator creates an empty 
state, whereas pseudo-fermion operator fj creates singly 
occupied state with an electron in the i-th dot. To elim- 
inate non-physical states, the following constraint has to 
be imposed on the new quasi-particles, 



Q = E///' + bt& = 1 - 



(6) 



Above constraint prevents double occupancy of the dots; 
dots are empty or singly occupied. 

In the next step the Hamiltonian |T]) of the system is 
replaced by an effective Hamiltonian, expressed in terms 
of the auxiliary boson b and pseudo-fermion /j operators 



as. 



H = 



E £ w c U c w 

k/3 



k L k 



E^/M + ^/i^ + H.c.) 

i 

H.c.)+A(X;4/< + 6 + 6-l J ■(?) 



k/3 i 



To avoid double occupancy of the DQD system the con- 
striction condition (Eq. (O) has been incorporated in 
Hamiltonian ||7J) by introducing the term with the La- 
grange multiplier A. 
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However, after such transformation our Hamiltonian 
is still rather complex and hard to solve. To get rid of 
this problem we apply mean field (MF) approximation 
in which the boson field b is replaced by a real and in- 
dependent of time c number, b(t) — > (b(t)) = b. This 
approximation neglects fluctuations around the average 
value (b(t)) of the slave boson operator, but is sufficient 
to describe correctly those leading to the Kondo effect. 
It also restricts our considerations to the low bias regime 
(eV< M)- 

With the following definitions of the renormalized pa- 
rameters: t = tb 2 , = Vflji and £, = £, + A, one can 
rewrite the effective MF Hamiltonian in the form, 

H MF = ]T e^cl^p + Y, Ziftfi + (tflh + kc.) 

k/3 i 

+EE(^U+ h - c -)+ A ( ? - 1 )- ^ 

k/3 i 

The unknown parameters b and A have to be found self- 
consistently with the help of the following equations; 

88 -*E/ = 1. (9) 

-<£/ ^{s-mmh)i + \b 2 = ^ (io) 

i 

where {{fi\fj))f is the Fourier transform of the lesser 
Green function defined as Gf ja (t,t') = ({fi{t)\f](t')))< = 
i(fj(t')fi(t)}. The above equations have been obtained 
from the constraints imposed on the slave boson, Eq. ([6]) , 
and from the equation of motion for the slave boson op- 
erator. The lesser Green functions ((fi\fj))f as well as 
retarded {{fi\f})) r (which also is needed in further cal- 
culations) have been determined from the corresponding 
equation of motion. 

To get insight into system's electronic properties we 
calculate local density of states (LDOS) and transmission 
function. The local density of states for the i-th dot is 
defined as; 

Di = --X>[G<i(e)]. (11) 

a 

whereas transmission through the system is expressed in 
the form; 

T{e) = Tr[G a f i? G r f L ], (12) 

In above formula T 13 stands for coupling matrix to the 
/3-th lead with renormalized parameters f£ = b 2 T^, and 
G r (G°) denotes Fourier transforms of the retarded (ad- 
vanced) Green functions of the dots. It is worth noting 
that transmission probability T(e) is directly related with 



current and linear conductance. In the zero temperature 
limit T — > these quantities are given by formulas 

e r v/2 

J=j deT(e), (13) 

n J-eV/2 

and 

^0=11^^ = ^ = 0). (14) 

III. NUMERICAL RESULTS 

A. Dicke effect 

In the following numerical calculations we assume 
equal dot energy levels, Ei = £q (for i = 1, 2) (eq is mea- 
sured from the Fermi level of the leads in equilibrium, 
(J'L = M-R = 0)- Moreover, we set the bare level of the 
dots at Eq = — 3r, and the bandwidth is assumed to be 
2D = 120r. All the energy quantities are expressed in 
the units of T. The parameters qi, and qr are assumed 
to be equal = qn — q if not stated otherwise. Tak- 
ing into account the above parameters, the Kondo tem- 
perature Tk for the symmetric couplings (a = 1) and 
disregarding both direct and indirect couplings, (t = 0, 
q = 0) is estimated to be Tk ~ 10~ 2 r. In this sec- 
tion we consider the situation when dots are not directly 
coupled, which corresponds to the case with vanishing 
interdot tunnel coupling parameter t = 0. On the other 
side, one should remember that dots are still coupled 
(indirectly) through the leads what is reflected in finite 
values of the off-diagonal coupling matrix elements, i.e., 
q =/= 0. In this situation bonding and antibonding levels, 
created due to indirect couplings, coincide. In Fig. [2] we 
show local density of states for QD1 and QD2 for large 
off-diagonal matrix couplings (q = 0.99) and for differ- 
ent values of the asymmetry parameter a. Firstly, it is 
clearly presented that the broad and narrow Kondo peaks 
in the LDOS are superimposed at energy e = and the 
LDOS displays behavior typical for the Dicke effect. By 
analogy to the original Dicke phenomenon, one may as- 
sociate the narrow (broad) central peak in LDOS with a 
subradiant (superradiant) state. The subradiant (super- 
radiant) state corresponds to longlived (shortlived) state. 
For symmetric coupling a = 1 the LDOS for QD1 and 
QD2 are the same (see Fig. HJa) but when asymmetry 
appears in couplings (a^ 1) this ceases to be true and 
the LDOS for both dots have different line shape. As a 
drops down the widths of the Kondo peaks for both dots 
also diminish. With decreasing a the broad part of the 
LDOS for the quantum dot strongly coupled to the leads 
(QD2) becomes more pronounced, whereas the narrow 
one becomes narrower. The two peaks corresponding to 
subradiant and superradiant state are well distinguish- 
able. As a result, the Dicke effect in LDOS for QD2 is 
more distinct. The Dicke effect is also noticed in the 
transmission shown in Fig.[3J This situation is analogous 
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FIG. 2: Local density of states for the quantum dots QD1 
and QD2 obtained for indicated values of a and for q = 0.99. 
Here, direct hopping between the dots is not allowed, because 
t = 0. 
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FIG. 4: Local density of states for the both dots (a) Friedel 
phase (b) and the transmission (c) calculated for indicated 
values of the q and for a — 1 and for t = 0. The well- 
defined Dicke line shape is only preserved for indirect coupling 
parameter q close to 1. 



FIG. 3: The transmission probability calculated for indicated 
values of a and for q = 0.99 and for t = 0. The transmission 
probability has well-defined Dicke line shape. 



to that reported in the case of the spin Kondo effect in 
parallel double dot system^ The Dicke peak appears 
in the transmission because the phases of the transmis- 
sion amplitudes for bonding and antibonding channels 
are equal at zero energy. Thus, the two contributions 
adds constructively leading to the maximum transmis- 
sion at e = 0. However, now there is no dip structure 
at zero energy for q = 1 (originating from the complete 
destructive interference), which will be explained further. 



Moreover, the effect is preserved for various values of the 
asymmetry parameter a. It is worth noting that the lin- 
ear conductance (see Eq. (fl"4)) ) reaches unitary limit (two 
quanta of e 2 /h). When a is reduced the effective Kondo 
temperature also decreases what can be seen looking at 
the energy scale in Figure [5J The origin of this effect 
comes from the fact that when a decreases one of the 
dots becomes detached from the electrodes. Then, the 
rate of the higher order tunneling events (which leads to 
the Kondo anomaly-see explanation in the Introduction) 
also diminishes and finally for a = 0, when one of the dot 
is totally disconnected from the leads, there is no possi- 
bility for such events and no Kondo effect is expected. 
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As we mentioned in Section III Al the off-diagonal ma- 
trix elements T^ 2 may be significantly reduced, so it is 
desired to analyze this case. In Fig. (E|) wc plotted the 
LDOS, Friedel phase and the transmission for various val- 
ues of the parameter q which is directly related with am- 
plitude of the off-diagonal matrix elements. The Friedel 
phase is related to the LDOS by the following equation: 
dOp/de = TvD(e) with D{e) being relevant density of 
states. One can notice that Dickc effect in LDOS can be 
found only when off-diagonal matrix elements are large, 
i.e., q close to 1. With decreasing q the Dickc line shape 
is transformed in usual Lorenzian line. Similar behavior 
is observed in the transmission (see Fig. Etc)). In inset 
of Fig. Etc)) we show that the narrow part of the peak 
has also well defined shape as the broad one. Moreover, 
when q is closer and closer to 1 then the narrow peak 
becomes more and more narrower. This is reflected in 
abrupt (but continuous as long as q ^ 1 ) change by 7r 
of Friedel phase around e — for q close to 1. It is also 
found that the Dicke effect also disappears in transmis- 
sion when q is maximal (q = 1). When q = 1 the trans- 
mission probability has Lorenzian shape and is described 
by formula, 



0.012 



2 I (l + a) 2 r 2 + (£ -e) 2 



(15) 



where f = b 2 T/2. This equation clearly shows that no 
Dicke effect should be expected for q = 1. This result 
resembles that obtained for a noninteracting system^! 
However, in the Kondo regime the situation is much 
more complex. To show this, let us first consider the 
case with q = 1 and symmetric couplings, a = 1. It is 
well known that for symmetric (a = 1) noninteracting 
(U = 0) system as q tends to 1 one of the peaks be- 
comes progressively narrowed and finally for q = 1 a BIC 
emerge s 34 ' 54 . As a result the transmission reveals simple 
lorenzian lineshape. One may naively believe that similar 
situation occurs in the Kondo regime. However, this can 
not be true because when the indirect coupling strength 
is equal to the dot-lead coupling (T 12 = Ln = L22), 
the indirect tunneling processes completely destroy co- 
herent higher-order tunclling events leading to complete 
suppression of the Kondo resonance^ 3 - One can also look 
at this from the another point of view and explain this 
as follows: As the antibonding state becomes a BIC, it is 
totally decoupled from the leads and there is no possibil- 
ity for an electron to exchange between the two molecular 
states. Thus, no Kondo effect appears as the correspond- 
ing Kondo temperature is equal to zero. To support these 
predictions we calculated the Kondo temperature for ar- 
bitrary value of the parameter q. The Kondo tempera- 
ture for the symmetric case, a = 1, acquires the follow- 
ing form Tk = \Jsq + f 2 (1 — q 2 ) 2 ■ In the deep Kondo 
regime the renormalized parameter £0 is equal to zero, 
which is shown in the Appendix, and thus the above for- 
mula clearly show vanishing of the corresponding Kondo 
temperature as q tends to 1. In Fig. [5] the Kondo tem- 
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FIG. 5: The Kondo temperature calculated as a function of 
parameter q and renormalized width T+ (T-) coresponding 
to the bonding (antybonding) level calculated for a = 1 and 
for t = 0. 
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FIG. 6: Graphical illustration of the Equation llGI The trans- 
mission probability calculated for indicated values of the a 
and for q — 0. 



perature is displayed as a function of parameter q. For 
comparison we also plot the q dependence of the renor- 
malized widths of the bonding and antibonding levels, 
and T 2 ., respectively. In the case q^ — qn — q, the widths 
acquire the following form r& = b 2 Ti,, T a = b T a with 
T h , a = §(T U + T 22 ) ± q^YuT^ 2 and r« = F£ + rg for 
i = 1,2. The characteristic widths for both distinct chan- 
nels behave in different way with varying the strength of 
the off-diagonal tunneling processes. One can notice that 
the renormalized width for the bonding channel changes 
non-monotonically but rather slowly, whereas the one for 
the antibonding channel drops monotonically to zero as 
q reaches maximum value. This behavior is consistent 
with the predictions given in Ref. [55| with one exception. 
The vanishing of T 2 when q is maximal is responsible for 
the suppression of the Dicke effect. In contrast to Ref. l55l 
the width Tg drops to zero for q = 1 due to vanishing of 
the boson field. We also emphasize that for q <G (0,1), 
the a dependances of the renormalized widths for both 
channels are different. For large values of the q the width 
T a changes little with a. 

However, the Kondo effect is also destroyed in more 



FIG. 7: The Kondo temperature calculated as a function of 
the asymmetry parameter a for q — and for t = 0. 
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general case, i.e., for q = 1 and arbitrary a. This can be 
understood when one notice that transmission (|T5|) de- 
pends on renormalizcd parameter b which vanishes in this 
case. In the SBMF formalism this is manifested as lack 
of solutions of the self-consistent equations ((9j) and (fT0|) 
(which then form a contradictory system of equations). 

On the other side, for q = the following formula 
describes the transmission probability, 



T{s) = \ 



r 2 + (e - e) 2 a 2 r 2 + (e - ef 



(16) 



Figure |6] illustrates graphically Eq. (fPoj) for indicated val- 
ues of the asymmetry parameter a. This plot clearly 
shows, as explained earlier, that the width of the Kondo 
peak decreases as a is reduced. This case resembles the 
spin Kondo effect in single-level quantum dot coupled to 
ferromagnetic leads^. Then the asymmetry parameter a 
can be assigned with lead's (pseudo)polarization in the 
following way p = (1 — q)/(1 + a)££. Increasing pseu- 
dopolarization (decreasing a) the Kondo effect is sup- 
pressed which is in agreement with Rcfs. [5611571 In the 
presence of asymmetry in couplings (a < 1) the Kondo 
temperature should be defined by geometric mean as: 
Tk = VT if 2 = T\/a (remembering that So — > in deep 
Kondo regime). Fig. [7] shows a dependence of the Kondo 
temperature for case of q = 0. In the present case, q = 0, 
electrons may travel only directly through dot QD1 or 
QD2 thus the reduction of the Kondo temperature with 
decreasing a clearly reflects the suppression of Kondo 
fluctuations in the DQD (due to reduction of the charge 
exchange between the dots) as one of the dot becomes 
detached from the leads and Kondo effect diminishes . 

Similar dependance on parameter q can be noticed in 
the differential conductance which is displayed in Fig. [5] 
The differential conductance is symmetric in respect to 
the zero bias point, thus we plot only its bias dependance 
for nonnegative bias voltage. For q close to unity the 
differential conductance acquires Dicke line shape. At 
zero bias the differential conductance approaches unitary 
limit (2e 2 /h) for all q < 1. For q close to 1 the differential 
conductance drops very fast to the half of the zero bias 
value with increasing the bias voltage and next decreases 



FIG. 8: Differential conductance as a function of bias voltage 
calculated for indicated values of parameter q (qL = qa, — q) 
and for a = 1, t = 0. 
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FIG. 9: The transmission coefficient calculated for indicated 
values of the parameter q {qL = —qn = q) and for a = 1. The 
inset show corresponding LDOS which is independent on the 
value of q. 



slowly with further increasing of the bias voltage. Such 
a sudden drop of the differential conductance is absent 
for smaller values of the parameter q. This feature in the 
differential conductance is related to Dickc resonance. 

As we mentioned in Sec. fl] the off-diagonal elements 
of the coupling matrix can have opposite signs. Specifi- 
cally, we examine the case qh = 1 and qn = —q. At the 
beginning we assume symmetric couplings, i.e., a = 1. 
Then for maximal value of the parameter q (q = 1) 
each DQD's molecular state couples to different reser- 
voirs and as there is no connection between them, the 
current doesn't flow. This effect originates from the 
totally destructive quantum interference which leads to 
vanishing of the transmission even if the LDOS of each 
dots is finite. It is wort noting that both LDOS and 
the corresponding Kondo temperature do not depend on 
the parameter q, as the contributions from nondiago- 
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nal processes, rcfcrcing to the left L and right R lead, 
cancel each others. This implies that the presented ef- 
fect originates fully from the quantum interference. One 
can notice that decreasing the value of the parameter q, 
the transmission can be recovered as is shown in Fig. [9] 
Finally, for q = the maximal value of the transmis- 
sion is restored. This is because as the value of the 
off-diagonal matrix elements are decreased the destruc- 
tive interference becomes totally suppressed, and thus, 
for zero value of the nondiagonal couplings the zero bias 
conductance is fully restored. This can be shown more 
formally by performing the transformation of the dots 
operators to the bonding-like (dj, = l/\/2(di + efe)) and 
antibonding-like (d a = \/y/2{d\ — cfo)) state. The corre- 
sponding couplings to the left lead acquires form a = 

5 (f 11 + ^22) ± Q'v/T uf 22' whereas couplings to the right 

electrode is given by f*. = |(ff x + f£) T qyfr&f*. 
Now, it is clear that for a = 1 and q = 1 the bond- 
ing state is coupled to (decoupled from) left (right) lead, 
whereas antibonding state is coupled to (decoupled from) 
right (left) electrode. As the parameter q becomes less 
than 1, the decoupled states from given leads for q = 1 
start to bound with them. As a result the transmission 
grows with decreasing the value of the parameter q. 
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FIG. 10: The transmission coefficient calculated for indicated 
values of the parameter a and for qL = —qn = 1. 



Another interesting quantum interference effect can be 
found changing the asymmetry in couplings of two dots 
to the leads. At the beginning we keep q equal to 1 
and change the parameter a in the interval (0,1). As 
before the transmission becomes recovered as the asym- 
metry increases [sec Fig.[TU]. However, in comparison to 
the previous case, a new feature emerges in the trans- 
mission. More specifically, the dip structure appears in 
the vicinity of the zero energy. At e = the transmis- 
sion drops to zero which results in vanishing of the linear 
conductance. Thus, the zero bias Kondo anomaly is to- 
tally suppressed for any value of the asymmetry param- 
eter a. To verify this effect experimentally it is desired 
to measure the differential conductance as a function of 
bias voltage. In Fig. [TT]wc displayed bias voltage depen- 
dence of the differential conductance for indicated values 
of the asymmetry parameter a. These results show that 
applying finite bias voltage, the differential conductance 
(which is zero at zero bias for all a) is restored. The 
differential conductance becomes the most pronounced 
when one of the dot is almost decoupled from the leads, 
i.e., as a — > 0, the differential conductance tend to uni- 
tary limit. However, with increasing asymmetry in the 
couplings (decreasing a) the range of finite values of the 
differential conductance shrinks due to decreasing of the 
Kondo temperature. It is worth noting that the corre- 
sponding Friedel phase is a continuous function of energy 
- it does not suffer discontinuity at e = because no BIC 
appears. However, the transmission at e = vanishes, 
thus the phase of transmission amplitude should be dis- 
continuous at e = 0. 
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FIG. 11: Differential conductance as a function of bias voltage 
calculated for indicated values of the parameter a and for 
q L = —qr = 1, t = 0. 



B. Fano effect 

Here, we show the results obtained for nonzero intcrdot 
hopping parameter (t ^ 0). In present situation the in- 
tcrdot hopping term lifts the degeneracy and the bonding 
and the antibonding levels are split away. As a result, the 
density of states of the DQD system coupled to the leads 
consists of two Kondo peaks; a broad peak centered at 
the bonding state and a narrow one corresponding to the 
antibonding state. This is clearly showed in Fig. [T^a)) 
where the LDOS for both dots is plotted for maximal off- 
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diagonal matrix elements (q = 1). It is very interesting 
that the narrow peaks in LDOS for QD1 and QD2 have 
opposite symmetries. The LDOS for QD1 reaches zero 
for negative energy, whereas the one for QD2 comes to 
zero for positive energy. Usually, when the asymmetry 
in coupling of two dots to the leads is reduced the width 
of bonding (antibonding) state increases (decreases). Fi- 
nally, for a — 1 (full symmetric system) the antibonding 
states is decoupled from the leads and acquires 5-Dirac 
shape (it becomes a BIC), whereas the bonding state ac- 
quires the width of 2T. On the other side, when the 
asymmetry in coupling of two dots to the leads is in- 
creased then the width of bonding (antibonding) state 
decreases (increases) and for a = the two resonances 
acquire the same width. However, in current problem the 
situation is more complex because the level's widths are 
renormalized by factor b which has to be obtained self- 
consistent ly for each a. Now it is only true that relative 
width of the bonding state to the width of the antibond- 
ing state grows as a increases, and are the same for a = 0. 
However, due to widths renormalization both peaks in 
LDOS may have smaller widths for small a with com- 
parison to the width of the narrow peak for large a (but 
not a ~ 1) (for instance for a = 0.15 both resonances are 
merged into narrow resonance for a = 0.8). Similar be- 
havior is observed in the LDOS for dots QD1 and QD2. 
Another feature is that for smaller values of the a only 
LDOS for QD2 may reach zero value at some energy. 
The above behavior of LDOS for a £ (0, 1) results in Fano 
antiresonance in the transmission probability. The cor- 
responding transmission probability calculated for q = 1 
and a = 0.8 is displayed in Fig. |T2fb) . One can notice 
that the transmission reveals the antiresonance charac- 
ter with a characteristic Fano line shape. In turn, the 
transmission function corresponding to the bonding state 
is relatively broad and (roughly) Lorcnzian. In inset of 
Fig. HHb) the T(e) is plotted for greater asymmetry in 
coupling of two dots to the leads (a = 0.15). It is no- 
ticed that the antiresonance behavior is still preserved 
but now the widths' ratio of the narrow resonance with 
respect to the broad one is much larger in comparison to 
the previous case (i.e., when a = 0.8). With decreasing 
a the widths of both peaks also decrease which can be 
seen looking at the energy scales for both cases. This 
means that the Kondo temperature is also lowered as a 
decreases (as was shown in the previous section). For 
q = 1 the transmission probability has poles at, 



2^ 

£q — — 1. 



1 



Q 



(17) 



It is worth to mention that at this energy the phase of 
the transmission amplitude suffers discontinuity. This 
is no more true for q < 1, because then the transmis- 
sion probability has complex poles. In turn, for q < 1 
the transmission probability does not reach zero and for 
specific value of q antiresonance behavior is less visible. 
Even a small shift in q change the transmission proba- 
bility within the antibonding level considerably, which is 
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FIG. 12: Local density of states for the both dots (a) and the 
transmission probability (b) calculated for q = 1, a = 0.8 and 
for t — 0.8. The existence of the Fano-Kondo effect is well 
visible. The inset shows the transmission probability obtained 
for q — 1 and for a — 0.15. 
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FIG. 13: The transmission probability calculated for a — 0.8, 
t = 0.8, and for q = 0.99. Even a small shift in q leads to the 
considerably changes in the transmission probability within 
the antibonding level considerably (compare with Fig. I12p . 



showed in Fig. 1131 

At this point we should remark that the slave boson 
mean field approach does not take into account the level 
renormalization arising due to coupling of dots to the 
leads. Such a renormalization should lead to splitting 
of the zero bias anomaly for specific cases^ The level 
splitting can occurs due to asymmetry in coupling of the 
dots to the leads, i.e., for a < 1 as has been shown in 
Ref . l56l by means of the scaling technique. However, such 
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a splitting can also be induced by indirect coupling of 
the dots. One can show that the splitting is proportional 
to the strength of the indirect coupling. Moreover, the 
level renormalization of a given state is proportional to 
the coupling strength of the other level. This implies 
that for q close to 1 the bonding-like level will be only 
weakly renormalized, whereas the antibonding-like level 
will experience strong renormalization. 

Thus, the transmission described in Sec. IIII Al should 
resembles that from Fig. rj5] but with the broad peak 
pinned close to e = and the narrow maximum shifted 
up in energy for relatively large value of parameter q. 
One can show this more formally including corrections in 
molecular-like levels due to mentioned renormalization. 

To correct the drawback of the SBMF method one 
can by hand introduce the mentioned renormalization of 
the levels as follows: Sb = £q + <5e& and e a = £o + Se a 
(with t = for the sake of simplicity). Here, <fej are 
the corrections due to indirect coupling of the dots lev- 
els. These corrections should be determined using rel- 
evant technique, as for instance earlier mentioned scal- 
ing procedure. However, as the SBMF technique fails 
for nondegenerate states, one can not introduce by hand 
the renormalized levels into self-consistent equations of 
the form @ and (JTUJ) derived within bonding and anti- 
bonding states basis. Thus, we give here only rough es- 
timation and some predictions implying from the levels 
scaling. Such approach should deliver qualitatively good 
insight into Kondo peak splitting phenomenon, however 
to obtain quantitatively consistent results more reliable 
technique should be applied for the considered problem. 
For the sake of simplicity we analyze only the symmetric 
case (a = 1). Introducing the level renormalization the 
transmission acquires the following form: 

= f 2 (l-g) 2 f 2 (l + g) 2 

(e - £a) 2 + f 2 (1 - q) 2 (e - e b ) 2 + f 2 (l + qf 

(18) 

with renormalized bonding (b) and antibonding (a) levels 
of the following assumed forrn^: 

Ei = e + TjA (19) 

Here, A££ stands for some function which in general de- 
pends on the system's parameters (like bandwidth, dots 
energy levels, couplings). For the sake of simplicity we 
assume A to be a constant number and being maximum 
value of function A(q) obtained from the scaling pro- 
cedure, i.e. when q = 1. Equation (|18[) clearly shows 
that Kondo peak becomes split due to the level renor- 
malization originating from indirect tunneling between 
the dots. In Fig. Q3] we display expected lineshape in 
the transmission for relatively large value of parameter 
q. The splitting in the transmission should decrease with 
decreasing the value of parameter q and beyond a certain 
value of q the splitting ceases to be visible. For q = 
(and a = 1) no splitting occurs. 



q=0.99 
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FIG. 14: The transmission coefficient including the level 
renormalization calculated for a = 1 and for qL ~ qn~ 0.99. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have investigated electronic proper- 
ties of double quantum dots coupled to external leads. 
The dots have been coupled both via hopping term and 
Coulomb interaction. Moreover, we have considered also 
the effects of indirect tunneling between the dots through 
the leads. Employing the slave-boson mean field ap- 
proach the local density of states for both dots, the 
Friedcl phase and the transmission in Kondo regime have 
been calculated. Moreover, to be more familiar with ex- 
periment we have calculated the corresponding differen- 
tial conductance. 

We have shown that for some set of parameters the 
Dicke-and Fano-like resonances may appear in the con- 
sidered system. More specifically, it has been noticed 
that for zero intcrdot hopping parameter t the LDOS of 
each QD consist of broad and narrow Kondo peaks that 
are superimposed. As in original Dicke effect one may as- 
sociate narrow (broad) central peak in LDOS with a sub- 
radiant (supcrradiant) mode. Moreover, Dicke line shape 
has been found in the transmission and in the differential 
conductance. We have observed that this effect is very 
sensitive to the strength of the off-diagonal matrix ele- 
ments; with reducing the value of the off-diagonal matrix 
elements the Dicke line both in LDOS and the transmis- 
sion is transformed into usual Lorenzian line. It has been 
found that when the intcrdot tunneling is allowed the 
transmission probability may reveal the antiresonance 
behavior with a characteristic Fano line shape. Moreover, 
we have noticed that the line shape of the antiresonance 
is also very susceptible to the change of the value of the 
off-diagonal matrix elements. 

We have also calculated the Kondo temperature, and 
also have shown that the latter becomes suppressed with 
increasing asymmetry in the dot-lead couplings when 
there is no indirect coupling. Moreover, when the in- 
direct coupling is turned on, the characteristic widths for 
both distinct channels behave in different way with vary- 
ing strength of the off-diagonal tunneling processes. We 
found also that the corresponding Kondo temperature is 
totally suppressed for maximal value of the indirect cou- 
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pling and no Kondo effect occurs. Moreover, we have 
also included level renormalization effects due to indirect 
coupling phenomenon, which leads to the splitting of the 
Kondo peak. 
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Appendix A: Proof of e — > in deep Kondo regime 

Here, we show that in deep Kondo regime the renor- 
malized parameter £o is equal to zero and, thus, the 
Kondo temperature strictly correspond to to the renor- 
malized width T of the Abricolov-Suhl resonance. This 
can be shown analytically by integrating self-consistent 
equations for slave-boson parameters, b, A, written in 
representation of bonding and antibonding states. 

In the basis of the bonding and antibonding states the 
Hamiltonian of the system becomes diagonal for e\ — 
£2 = £o an d acquires the following form: 



H = H c + Ejdldj + Un b n a 



i—b,a 

+ EE(^+ Ec -)> (ai) 

k/3 i— e,o 

with Eb = £o + t, e a = £o — t. In further considerations 
we assume t = 0, thus, £/, = s a = Eq. 

In the mean field slave boson representation the Hamil- 
tonian (|Alj) acquires the following form: 



H 



MF 



EEft/. + H-c). (A2) 

k/3 i—b,a 



The self-consistent equations determining the un- 
known parameters b and A have the form of Eqs. ([9]) 
and (fTU]) with lesser Green's function ((fi\fj))f defined 
in the basis of the bonding and antibonding states. The 
lesser bonding and antibonding Green's function has the 
following form: 



(e - £ ) 2 + (r,/2) 



(A3) 



At T — K f(x) = 6{—x) and the integration of the 
self-consistent equations is straightforward^ leading to 
the following expressions: 



l-b 2 




Xb 2 = 0. (A4) 



i—b,a 



In the deep Kondo regime we can approximate: 1 — b w 1 
and A w — Eq and then Eqs. (|A4|) simplify to the following 
equations: 




27T£ . (A5) 



Here, we consider only the symmetric case, a = 1, how- 
ever, extension to arbitrary a is straightforward. For 
a = 1 and qL = qn = q, the couplings to the bond- 
ing and antibonding states acquires the following form, 
Fb.a = (1 ± q)T. Combining two equations (|A5[) one ar- 
rives with equation: 



(e +if b /2)(e + it a /2) 



E + (f fc /2) 2 

^0 + (f o /2) s 



D exp 



£o 

r 



2k 



(A6) 



which real and imaginary parts satisfy the following equalities: 
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2 /'27T£ 

-D exp 



= 0. 



(A7) 
(A8) 



Equation (|A8|) is satisfied if eo = or 6 2 = 0. However, It is worth rioting that the above equation does not de- 
thc later solution leads to nonphysical value for io, thus, tcrmine the value of b 2 for q = 1 when the SBMF method 
the only solution must be io = 0. This result clearly fails, 
shows that in deep Kondo regime slave-boson parameter 
eo vanishes, Q.E.D. 

Introducing the solution for io into Eq. (|A7j) one finds, 

V=™ {l - q) Z ™A—)- (A9) 
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